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We extend the eV seesaw scenario to four lepton generations. The LSND anomaly is taken as the 
right-handed seesaw scale, i.e. tur ~ eV. The fourth generation then gives a heavy pseudo-Dirac 
neutrino which largely decouples from other generations, and is relatively stable. One effectively 
has a 3 -I- 3 solution to the LSND anomaly, where we illustrate with numerical solutions. Our study 
seems to indicate that the third mixing angle sin^ 6^13 may be less than 0 . 01 . 
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I. INTRODUCTION 

Neutrino oscillation measurements have become more 
and more accurate. The latest combined results [1] 
read Atuq = (7.9 ± 0.7) x 10“® eV^and = 

(2.4to;6) X 10“^ eV^, for the mass squared differences of 
the large mixing angle solution (LMA) to the solar neu¬ 
trino problem [2, 3], and for atmospheric neutrinos [4-6], 
respectively. Both mass differences are sub-eV, but the 
neutrino mass scale is not yet certain. 

Neutrino data also hint at the possibility of more 
than three massive, mostly active neutrinos. The Liq¬ 
uid Scintillator Neutrino Detector (LSND) Collabora¬ 
tion [7] has reported evidence for a flux 30 meters 
away from a source of produced in 7r+ —> 
with subsequent /i'*' -\- e'^ -\- Ve decay. The unex¬ 

pected flux can be explained if there is a small probability 
P^Pfj, Pe) = (0.26 ± 0.08)% for a neutrino produced as 
a Pfj, to be detected as a Pe [8]. This LSND anomaly has 
yet to be confirmed, and the MiniBooNE experiment at 
Fermilab [9] should very soon give a definitive confirma¬ 
tion or refutation. Numerous attempts to solve the LSND 
puzzle, however, have already been proposed [10]. It has 
been shown that oscillations with extra sterile neutrinos 
can fit the LSND anomaly [11]. But it has also been 
pointed out that extra sterile neutrinos could be in con¬ 
flict with Big Bang Nucleosynthesis (BBN) [12], as well 
as SN1987A supernova neutrino events [13]. 

In a recent work [14], it has been argued that, if the 
right-handed neutrino Majorana scale mu is of 0{eV), 
adequate fits to the LSND data can be obtained. This 
“eV seesaw” scenario runs against theoretical arguments 
in favor of a very large mn. To name a few such ar¬ 
guments: the canonical seesaw mechanism [15-17] with 
mji ~ 10^^ GeV can elegantly explain why neutrino 
masses are so small, even with lepton Yukawa couplings 
that are of order one; thermal leptogenesis [18] points to 
itir P 10^° GeV [19]. However, as stressed in Ref. [14], 
nothing experimental is really known about the magni¬ 
tude of ruR, except perhaps the LSND result, which is at 
eV scale. 

The purpose of this letter is to show that the eV seesaw 
proposed in Ref. [14] can be straightforwardly incorpo¬ 
rated in a four generation scenario. 


The Standard Model (SM) with a sequential fourth 
generation (SM4) is not ruled out by electroweak preci¬ 
sion measurements, if one allows the extra active neutrino 
to have mass close to 50 GeV [20, 21]. To avoid bounds 
from direct search at LEP II [22], mixing of the fourth 
heavy neutrino with the three light neutrinos should be 
small (^ 10“®). It is clear that in standard seesaw with 
rriR ~ 10^® GeV, an extra generation is hard to accom¬ 
modate (A different approach to predict a light sterile 
neutrino in the presence of a fourth generation is the so 
called “flipped seesaw” [23]). All four (mostly) active 
neutrinos will then be light, contradicting the invisible 
Z width which measures only three light neutrinos. But 
taking niR at scale ©(eV), one can now have a sufficiently 
heavy fourth neutrino. 

In the following we will show that, by taking mR ~ 
©(eV), the fourth neutrino is pseudo-Dirac and heavy. It 
will not affect the invisible Z width, and largely decouples 
from lower generations. Aside from three mostly active 
light neutrinos, three sterile neutrinos with mass ^ eV is 
predicted. A numerical analysis gives results consistent 
with the LSND as well as solar and atmospheric data. It 
seems that sin^ 6*13 cannot be large. 


II. PSEUDO-DIRAC FOURTH NEUTRINO 

Following Ref. [14] but allowing for a possible 4th gen¬ 
eration, the 8x8 neutrino mass matrix M is given by 

M = MD + AMR + dMD, (1) 

in a form suggestive of mass hierarchies. In the basis 
where the 4x4 Dirac mass matrix is diagonal, the dom¬ 
inant Dirac mass for the 4th generation arises from 

= mz, ° ^ , (2) 

where tod ~ 50 GeV, 0 and I 4 are 4x4 matrices with 
zero elements, except 1 in 44 element of Iq. The right- 
handed Majorana mass matrix is given by 

AMr = mzi ^ q ® ^ , (3) 


where mu ~ eV [14], and r is a 4 x 4 symmetric matrix 
with elements ~ 1. The third matrix is 


5Md = mu 



( 4 ) 


which is pinned more to the mu scale, with 

/ ei 0 0 0 \ 

_ I 0 £2 0 0 I 

0 0 £3 0 ’ 

Vo 0 00 / 


( 5 ) 


where £4 has been absorbed into mo- Clearly, with 
mR/mo = X ~ 10 “^° and ei considerably less than 1 , 
/S.Mr and 5Md can be treated as perturbations to Mr. 

We note that the neutrino mass matrix of Eq. (1) could 
arise from very small deviations from a democratic struc¬ 
ture for the Dirac contribution [24], and lepton number is 
assumed to be only slightly violated. The latter requires 
~ 0 [14], i.e. symmetry is enhanced in the limit of 
AMr 0. Having smaller elements in 6Mr compared 
to mR is purely phenomenological [14]. 

The dominant Mr has six null eigenvalues, plus 
two eigenvalues ±m_D. For the eigenvectors associ¬ 
ated with zero eigenvalues, one can choose = 

( 0,..., 0, 1 , 0 ..., 0 ) with 1 in the position for i = 

1,2,3 and i = 5,6,7. For the eigenvalues —mR 
and mR, the corresponding eigenvectors are e 4 °g = 
( 0 , 0 , 0 , =f 1 /-\/ 2 , 0 , 0 , 0 , l/-\/ 2 ). At zeroth order in x, the 
states £ 4 °^ and e® combine into a pure Dirac state of 
mass mR. When linear corrections in x are considered, 
the perturbed states £4 and £§ with Ci = £^°^ -|-x£^^^ have 
masses which differ by 0 {x), and they now correspond 
to a pseudo-Dirac neutrino with mass ^ mR, which we 
denote as N (the charged partner is denoted E, with 
ms > 100 GeV [25]). 

For the six null eigenvalues ot Mr, one can apply per¬ 
turbation theory with degeneracies. A linear combina¬ 
tion of the degenerate unperturbed states £.°^ diagonal¬ 
izes the AMr -\- SMr perturbation. That is, by diago¬ 
nalizing the 6 x 6 perturbation matrix 

0 0 0 ei 0 0 \ 

0 0 0 0 £2 0 

0 0 0 0 0 63 

€1 0 0 rii ri2 ri3 ’ v^/ 

0 62 0 r21 r22 ^23 

0 0 €3 rsi r32 r33 / 

one obtains the corrections of 0 {x) to the mass eigenval¬ 
ues, and the correct eigenstates at order zero. For the 


= mR 


effect of the 4th generation neutrino through the right- 
handed sector, Tij, one has linear corrections in x pro¬ 
portional to r44 to the eigenvalues imu for £4 and £3 
as already stated, and ©(x^) corrections to all the other 
eigenvalues. The big hierarchy between the matrix ele¬ 
ments M 48 , M 84 = mR and all the others allow the fourth 
generation to largely decouple from the other three. As 
stated in the Introduction, this is also required by direct 
search limits that demand very small mixings between N 
and the light neutrino flavors. 

Having reduced the problem to a 6 x 6 case, the analysis 
performed in [14] suggests that one could find a solution 
to the LSND puzzle. Our main goal here is to confirm 
the possibility of an existing solution, and to gain some 
insight on what could be a plausible scenario. 


III. 3+3 NEUTRINO MODEL 


We set to zero all phases for simplicity, since there are 
already too many parameters. We define U' to be the 
rotation matrix which diagonalizes M^^\ Having started 
in the basis where the Dirac neutrino mass matrix Mr + 
SMr is diagonal, one still has the freedom to perform a 
rotation U” in the left sector 


U" 


/ C 1 C 3 

I -S1C2 — C1S2S3 

I S1S2 — C1C2S3 

0 
0 
0 


S 1 C 3 

C1C2 — S1S2S3 
— C1S2 — S1C2S3 
0 
0 
0 


S3 0 0 0 \ 

S 2 C 3 000 
C 2 C 3 0 0 0 

0 10 0 
0 0 10 
0 0 0 1 / 


( 7 ) 


We assume no mixing between the fourth and first three 
generation charged leptons, as already discussed. For 
the right sector, a rotation will just change to rb, 
resulting in no change to our numerical analysis. 

The probability for a neutrino, produced with flavor a 
and energy E, to be detected as a neutrino of flavor fH 
after travelling a distance L is [25] 


^ — ^a (3 4 ^ ] Ua^jUj 3 jU(xiUpi sin Xji, (8) 

j>i 


where a = e, /i, r, Si with Si the sterile neutrino flavors, 
U = U"U', and Xji = 1.27 Am^^^ L/E with Am|j = — 

mf. Applying Eq. ( 8 ) [26] in the “3 active plus 3 sterile 
neutrino” (3 + 3) case, using the approximations X 21 = 
2^31 = X32 = 0 and Xii = Xi 2 = Xi^ for i = 4,5,6, one 
obtains 


P{Ra vp) = Sap + 4[f7^4(f7^4 - Sap) sin^ X41 + Ul^{Up^ - Sap) sin^ X51 + UIq{Upq - Sap) sin^ xei 

+Ua 4 Up 4 Ua 5 Up 5 (sm^ X41 + sin^ X51 - sin^ X54) + Ua4Up4Ua&UpQ{sir? X 41 + sin^ X61 - sin^ X64) 
EUaiitUpQUa^UpQi^sin X51 + sin xqi sin X(55)], (9) 
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TABLE I: Best fit values, 2a and 3a intervals for three-flavor 
neutrino oscillation parameters from global data, including 
solar, atmospheric, reactor (KamLAND and CHOOZ) and 
accelerator (K2K) experiments, taken from Ref. [28]. 



BEST FIT 

2cr 

3cr 

Amii (10-® eV^) 

8.1 

7.5 - 8.7 

7.2 - 9.1 

Amii (10-® eV^) 

2.2 

1.7-2.9 

1 . 4 - 3.3 

sin^ 012 

0.30 

0.25 - 0.34 

0.23-0.38 

sin^ 023 

0.50 

0.38 - 0.64 

0.34 - 0.68 

sin^ 013 

0.0 

< 0.028 

< 0.047 


where orthogonality of U has been used. Expressions for 
the mixing angles are given by [27] 


tan^ 012 
sin^ 013 


\Uel\^ 


tan^ 023 = 


\Ur3? ’ 


|C/e3r 


( 10 ) 


To perform our numerical analysis, we build the 
by using the three-flavor neutrino oscillation parameters 
and sin^ 6 ij taken from Ref. [28], which is com¬ 
piled from global data including solar, atmospheric, re¬ 
actor (KamLAND and CHOOZ) and accelerator (K2K) 
experiments. These are given in Table 1. We also include 
the LSND result of P{v^i 7'e) = (0.26 ± 0.08)%, and 
require Am^i ~ leV^, for a total of seven inputs. As 
one clearly has too many parameters to make a proper 
fit, we just minimize the built with these quantities 
by letting the twelve parameters of the model vary. 

As an illustration, we find si = —0.57, S2 = 0.98, S3 = 
0.80, which give 


0.49 

-0.34 

0.80 

0 

0 

o\ 

-0.54 

0.60 

0.58 

0 

0 

0 

-0.69 

-0.72 

0.11 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 / 


( 11 ) 


and 




5eV 


( ° 

0 

0 

0.051 
0 
0 


0 0 

0 0 

0 0 

0 0 

0.004 0 

0 0.031 


0.051 0 0 \ 

0 0.004 0 

0 0 0.031 

1.195 0.861 1.038 
0.861 0.968 0.878 
1.038 0.878 1.264 / 


( 12 ) 

from which, together with Eq. (6), one can read the val¬ 
ues for the rest of the parameters coming from the mini¬ 
mization process. The eigenvalues of are 


(mi, m2, m3) = —(7 X 10 9 X 10 0.048) eV, 

(m4, ms, me) = (1, 1.2, 15) eV, (13) 


and the associated rotation matrix is 


/ -0.06 

-0.99 

-0.11 

0 

0 

0 

-0.38 

0.12 

-0.91 

0.01 

-0.06 

0.05 

0.90 

-0.02 

-0.37 

-0.17 

0.05 

0.12 

-0.19 

0 

0.09 

-0.75 

0.16 

0.60 

0.05 

-0.01 

0.07 

0.22 

-0.83 

0.49 

\ 0.01 

0 

0 

0.60 

0.52 

0.61 


(14) 


From Eq. (13) one gets Am^i = 8.1 x 10“^ eV^ and 
Am|i = 2.3 X 10“^ eV^, which of course is in good 
agreement with data. The requirements for the neutrino 
mass splitting applied in Ref. [26] for the 3-1-2 case, 
0.1 eV^ < Am4;^ < Amg^ < 100eW, are also satisfied. 
This guarantees that the approximations used to derive 
Eq. (9) are valid. 

From Eqs. (11) and (14), we obtain the full rotation 
matrix U = U"U', i.e. 


/ 

0.82 

-0.54 

-0.04 

-0.14 

0.06 

0.08 

\ 


0.33 

0.60 

-0.71 

-0.09 

-0.01 

0.10 


0.42 

0.59 

0.69 

-0.03 

0.05 

-0.03 



-0.19 

0 

0.09 

-0.75 

0.16 

0.60 


1 

0.05 

-0.01 

0.07 

0.22 

-0.84 

0.49 

/ 

0.01 

0 

0.01 

0.60 

0.52 

0.61 


Using Xji = 1.27Am|j(eV^)T(m)/if(MeV) with L/E ^ 
1, together with Eq. (9) for the fi ^ e case, Eqs. (13) 
and (15), one obtains P{v^ De) = 0.15%, which is 
within 2a from the LSND central value. 

in Eq. (12) can be viewed as deviating from 




ruR 


/oooooo\ 
0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 1 1 1 

0 0 0 1 1 1 

\ 0 0 0 1 1 1 / 


(16) 


which has five zero eigenvalues and one nonzero eigen¬ 
value equal to 3>mR ~ 15 eV, and diagonalized by 


U' = 


/I 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

n/2 

1 

0 

0 

0 

0 

1 

0 

1 

\0 

0 

0 

1 

x/3 

1 

v/3 

1 


Together with U" of Eq. (7), one has 


(17) 


C1C3 

S1C3 

S3 

0 

0 

0 

\ 

—S1C2 — C1S2S3 

C1C2 — 5152-53 

S2C3 

0 

0 

0 


SlS 2 — C1C2S3 

— ClS2 — S1C253 

C2C3 

0 

0 

0 


0 

0 

0 

1 

1 

0 


x/2 

x/2 


0 

0 

0 

_^ 

0 

1 


V 2 

\/2 

/ 

0 

0 

0 

1 

V3 

1 

V3 

1 


(18) 


Applying this to Eqs. (9) for the fi ^ e case, one sees 
that the transition probability P{i'p, —>■ Ve) would van¬ 
ish. Deviations from Eq. (16) as realized in Eq. (12) not 
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FIG. 1: Contour-plot of vs the mixing angles si and S 3 , 
with ei and Vij as in Eq. (12), and S 2 = 0.98 held fixed. The 
regions in different shades are only indicative, and should not 
be interpreted as the Ict, 2cr and 3a regions, as the rest of the 
parameters are fixed at the best fit values. 


P(z/e ^ i^e) = 0.89 and ^ = 0.93 in the LjE ~ 

1 approximation. Our calculated values for oscillation 
appearance and disappearance probabilities can now be 
compared with the ones obtained by using the best fit 
values for the rotation matrix elements in the (3-1-2) case 
as in Ref. [26], = 0.121, [/^4 = 0.204, t/es = 0.036 and 

U)i 5 = 0.224. In the LjE ~ 1 approximation these give 
P{v^i Ve) = 0.0021, P{ve Ve) = 0.95 and P{Vfi 
v^) = 0.84. We remark that, although a full analysis is 
needed to tell if our model is able to accommodate both 
NSBL and LSND data, our predictions seem to be not 
too far away from the results of Ref. [26]. 

From Eqs. (10) and (15), we find sin^ 012 = 0.30, 
sin^ 023 = 0.52, and 

sin^ 013 = 0.0018. (19) 

As expected, the values for sin^ 0i2 and sin^ 023 are in 
good agreement with data, but the as yet unmeasured 
sin^ 013 turns out to be rather small. 


only produces the needed neutrino mass spectrum, finite 
transition probability P{vfi —> Ve) is also achieved. 

In our analysis we did not take into account the null 
short-baseline experiments (NSBL). Let us compare our 
results with the ones obtained by doing a full analysis in 
Ref [26] for the (3-1-2) case. Using the rotation matrix 
elements t/e4 = —0.14, = —0.09, Ue 5 = 0.06, = 

—0.01, Uee = 0.08 and = 0.10 as one can read from 

Eq. (15), together with Eq. (9) respectively for /r ^ e, 
e ^ e and fi —>■ fi, we obtain P{iy^ Ve) = 0.0015, 


IV. IN SEARCH OF SIZABLE sin^ 0i3 

To investigate the possibility for a bigger value of 
sin^ 013, we restrict the to the four inputs of sin^ 6 ij 
and P{i>ii Ve)- The mass spectrum is not affected 
by the rotation of Eq. (7). With and given as in 
Eq. (12), we first fix si = —0.57 and perform a fit vs 
S2 and S3. We iterate with fixing S2 = 0.98 (S3 = 0.8) 
and minimize x^ vs Si and S3 (si and S2). We find for 
both cases of fixing si and S3 to the values found in pre¬ 
vious section, sin^ 023 is quite strongly dependent on S2, 



FIG. 2: P{Vfi —> Ve), sin^ 012, sin^ 023 and sin^ 0i3 vs Si and 
S3, corresponding to the lower right solution in Fig. 1, with 
ei and rij fixed as in Eq. (12), and S 2 = 0.98. 
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FIG. 3: Same as Fig. 2, but corresponding to upper left 
solution in Fig. 1. 









and the value around 0.98 is preferred. We thus illustrate 
with fixing S2 = 0.98. 

In Fig. 1 we show the contour plot of vs si, S3. The 
three different shaded regions should not be interpreted 
as the Icr, 2ct and 3 ct regions, since we have fixed the rest 
of the parameters to the best fit values. But they still 
give an indication of variations around the best fit region 
under the above assumptions. 

In Fig. 2 we plot the four quantities P{v^ f'e), 
sin^ 012, sin^ 023 and sin^ 0i3 vs si and S3, for the solution 
on the lower right of Fig. 1. The same is plotted in Fig. 3 
for the upper left solution. Again, and are fixed as 
in Eq. (12), and S2 is held fixed at 0.98. We see that 
Vg) can reach the one a region and sin^ 0i2 is 
well within range. However, to push sin^ 0i3 beyond 0.01, 
sin^ 023 seems to wander away from maximal mixing of 
0.5, and values at ~ 0.4 or 0.6 has to be tolerated. We 
note further that the sensitivity of sin^ 023 is to si, rather 
than S3. 

We conclude that sin^ 0i3 greater than 0.01 is possible, 
but seemingly not preferred. It is not clear whether this 
is an artefact of not being able to do a real fit. Note 
that we have not checked explicitly whether constraints 
from short baseline disappearance experiments are fully 
satisfied. 


V. DISCUSSION AND CONCLUSION 

It is tempting to consider whether mixing between the 
fourth and the first three light charged lepton generations 


could modify the situation with sin^0i3. But as already 
mentioned in the Introduction, one needs to satisfy both 
bounds from direct search at LEP II and electroweak 
precision measurements. We have pursued a numerical 
study, but find that, if we wish to keep the mixings suf¬ 
ficiently small so that the fourth active heavy neutrino 
will be semi-stable, no important change with respect to 
the no-mixing case is observed. We note that the heavy 
neutrino could be heavier than 50 GeV and still with sup¬ 
pressed mixing to lower generations, but then one would 
have to face electroweak precision constraints. We note 
in passing that semi-stable heavy neutrinos are still of in¬ 
terest [29] to dark matter search experiments, as a fourth 
heavy lepton was once a leading dark matter candidate. 

In summary, we have extended the eV seesaw scenario 
to four lepton generations. Taking the LSND scale as 
the right-handed seesaw scale mu ~ eV, one has a heavy 
pseudo-Dirac neutrino with mass rriAr ~ 50 GeV, which 
largely decouples from other generations, and is relatively 
stable. One effectively has a 3 -I- 3 solution to the LSND 
anomaly, where we illustrate with numerical solutions. 
As a possible outcome, our numerical study indicates 
that the third mixing angle, sin^0i3, seems to be less 
than 0.01. 
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